Nonlinear terahertz metamaterials via field-enhanced carrier dynamics in GaAs 
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We demonstrate nonlinear metamaterial split ring resonators (SRRs) on GaAs at terahertz fre- 
quencies. For SRRs on doped GaAs films, incident terahertz radiation with peak fields of ~20 - 
160 kV/cm drives intervalley scattering. This reduces the carrier mobility and enhances the SRR 
LC response due to a conductivity decrease in the doped thin film. Above ~160 kV/cm, electric 
field enhancement within the SRR gaps leads to efficient impact ionization, increasing the carrier 
density and the conductivity which, in turn, suppresses the SRR resonance. We demonstrate an 
increase of up to 10 orders of magnitude in the carrier density in the SRR gaps on semi-insulating 
GaAs substrate. Furthermore, we show that the effective permittivity can be swept from negative 
to positive values with increasing terahertz field strength in the impact ionization regime, enabling 
new possibilities for nonlinear metamaterials. 



Nonlinear metamaterials is a rapidly developing field 
of fundamental interest with significant technological im- 
plications spanning from microwave through the visible 
spectral ranges P~HI] ■ As with tunable and reconfigurable 
metamaterials [H [9] , the combination of the metamate- 
rial structure with the local environment is crucial. This 
is because significant nonlinearities result from local field 
enhancement within the active region of the subwave- 
length metamaterial elements which, in the case of split 
ring resonators, are the capacitive gaps. While the ac- 
tive volume of the enhanced gaps is small in comparison 
to the unit cell volume, the field enhancement can dom- 
inate volumetric effects leading to global nonlinearities 
enhanced by two to four orders of magnitude [5]. This, 
in turn, results in useful nonlinear effects at low incident 
fields. 

Advances in nonlinear metamaterials coincide with the 
development of high-field terahertz sources capable of 
generating electric fields sufficient to induce significant 
nonlinearities in conventional matter |10H14j . For ex- 
ample, in doped GaAs, highly nonlinear effects such as 
velocity saturation and impact ionization have been ob- 
served [HI [14] at peak electric fields of several hundred 
kV/cm. Further, with metamaterial THz field enhance- 
ment to MV/cm fields an insulator-metal phase transi- 
tion has been induced in vanadium dioxide, a prototypi- 
cal correlated electron material [7] . 

In this letter, we experimentally demonstrate a nonlin- 
ear response in metamaterial split ring resonators (SRRs) 
on n-type GaAs and semi-insulating (SI) GaAs at tera- 
hertz frequencies. The nonlinear response arises from 
THz electric field-induced carrier dynamics that increase 
or decrease the substrate conductivity upon which the 
SRR arrays are fabricated. This modifies the SRR elec- 
tromagnetic response as a function of field strength. For 
peak incident THz fields (E in ) from ^20-160 kV/cm, 



mobility saturation by intervalley scattering (IVS) dom- 
inates leading (for doped GaAs) to a conductivity de- 
crease and a corresponding increase in the metamaterial 
oscillator strength. In this regime, electric field enhance- 
ment within the SRR capacitive gaps does not signifi- 
cantly contribute to the nonlinearity. That is, IVS occurs 
uniformly across the film independent of the SRRs. For 
Ei n > 160 kV/cm, field enhancement in the gaps enables 
efficient carrier generation via impact ionization (IMI) 
resulting in a decrease in the SRR oscillator strength. It 
is possible to sweep the effective permittivity from neg- 
ative to positive values by increasing the incident field 
strength enabling new possibilities for nonlinear metama- 
terials such as an intensity dependent effective refractive 
index. We also show that SRRs on semi-insulating GaAs 
exhibit a nonlinear response arising from an increase of 
the in-gap conductivity by several orders of magnitude. 
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FIG. 1: (a). Image of the SRR array on n-type GaAs sub- 
strate. The period p is 50 /jm; the length of ring I is 36 [im; 
the gap g is 2.2 [im. (b) Field transmission spectrum of SRR 
array using THz time domain spectroscopy on doped GaAs 
at low fields (i.e. in the linear response regime). 



2 



To create nonlinear metamaterials, electric field reso- 
nant SRR arrays were fabricated by conventional micro- 
fabrication on GaAs. Fig. 1 (a) shows an SRR array on a 
1.8-^m-thick doped (f x 10 16 cm -3 ) GaAs film deposited 
on a SI-GaAs substrate using molecular-beam epitaxy. 
The SRRs consist of 200 nm thick Au with a 10 nm Cr 
adhesion layer. A narrow gap (g=2.2 /jm) is important as 
the field enhancement scales as 1 / g [15] . To characterize 
the linear electromagnetic (EM) response, electro-optic 
based terahertz time-domian (THz-TDS) spectroscopy 
was employed. To accurately measure the metamate- 
rial transmission spectrum, a bare n-type GaAs substrate 
was used as a reference. Dividing the Fourier transform 
of the sample data by that of the reference data yields 
the transmission amplitude as a function of frequency 
as shown in Figure [TJb). The LC resonance is ~ 0.75 
THz and is strongly damped, as expected, since carriers 
in the doped GaAs short the SRRs. The corresponding 
conductivity of doped GaAs is 7 (ficm) -1 . 

For the nonlinear transmission measurements, high- 
field THz experiments were performed using tilted-pulse- 
front THz generation in LiNbOs [TB] . We obtained a peak 
THz field strength of approximately 400 kV/cm, with a 
1 mm beam diameter at the focus. The field strength 
could be varied from 24 kV/cm to 400 kV/cm using a 
pair of wire grid polarizers. The detailed experimental 
setup has been described elsewhere [TUl E] ■ All experi- 
ments were performed at room temperature in a dry air 
environment at ambient pressure. 

Figure^ a) shows the experimentally measured trans- 
mission as a function of frequency for various values 
of Ei n < 160 kV/cm. At 24 kV/cm the transmission 
response is similar to the low-field (<1 kV/cm) result 
shown in Fig. [ljb) where a weak LC resonance is ob- 
served. As the incident field is increased to 160 kV/cm, 
the LC resonance becomes quite pronounced leading to 
a decrease in the resonant transmission (~ 0.73 THz) 
from 65% to about 30%. This suggests a decrease in 
the conductivity which, as discussed below, results from 
a decrease in the carrier mobility due to THz-induced 
I VS. The off-resonance transmission at low frequencies 
increases (e.g. at 0.4 THz from 80% to over 90%). This 
indicates a field-induced reduction in the conductivity 
of the entire n-type GaAs film similar to other tunable 
metamaterial responses [THJ, [Jj5] and consistent with ear- 
lier observation in bulk GaAs [TT]. Thus, for 24 < Ei n 
< 160 kV/cm the nonlinearity does not depend strongly 
on the in-gap field enhancement. 

However above ~ 160 kV/cm, the trend described 
above in the field-dependent transmission at the reso- 
nant frequency is exactly reversed as shown in Fig. 2(b): 
the transmission decreases with further increase in Ei n . 
This indicates a reversal of the underlying trend in con- 
ductivity as well, i.e. an increase in conductivity in the 
gaps as the field is increased above 160 kV/cm. Notably, 
the off-resonance transmission at lower frequencies (< 0.5 



THz) remains almost unchanged and implies a local con- 
ductivity increase solely in the capacitive gaps of the res- 
onators. We ascribe this response to THz field-induced 
impact ionization. 

To better understand how THz field-driven nonequi- 
librium carrier transport determines the nonlinear meta- 
material response, we consider, in turn, IVS and EMI. 
For 24 < E in < 160 kV/cm, electrons in the conduction 
band (r valley) are accelerated and acquire energy on 
the order of 1 eV, sufficient for scattering into satellite 
valleys. Since the mobility in the side valleys is approx- 
imately one order of magnitude smaller than in the T 
valley [2U1 HI]> IVS leads to a significant conductivity 
decrease in the n-type GaAs film. Independent of the 
SRR array, for doped substrates IVS results in photoin- 
duced transparency [TTJ [22l [23] and accounts for the off- 
resonance transmission increase that is observed in Fig. 
2(a). Of course, for the SRR arrays, this same process de- 
creases the in-gap conductivity, resulting in an enhanced 
LC resonance. 

IVS-induced enhancement of the metamaterial EM re- 
sponse was simulated using a commercial software, CST 
Microwave Studio. The n-type GaAs film was modeled 
utilizing a frequency dependent Drude model. We as- 
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FIG. 2: (a) & (b) Experimental electric field transmission 
as a function of frequency for various incident peak electric 
fields, (c) Simulated transmission of the SRR array on doped 
GaAs film for different electron mobilities at constant carrier 
density of 1 x 10 16 cm~ 3 . In these simulations the mobility 
of the carriers throughout the film is changed, /j, is the elec- 
tron mobility, (d) Simulated transmission on doped GaAs for 
various in-gap conductivities. In these simulations the mobil- 
ity is constant and the carrier density is modified according 
to IMI simulations. The inset in (d) shows the simulated 
on-resonance field enhancement highlighting that the largest 
enhancement occurs within the SRR gaps. 
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sume that the average mobility varies by changing the 
relative fraction of electrons in the T and L valleys and 
that at 160 kV/cm all of the conduction electrons are 
scattered to the L valley. Additionally, for fields < 160 
kV/cm, we neglect IMI. That is, the carrier concentra- 
tion is held constant at the doping value (1 x 10 16 cmT 3 ). 
The results of this mobility averaging procedure (applied 
to the entire film) determine the film conductivity used in 
the electromagnetic simulations. The results are shown 
in Fig. 2(c) and show good agreement with the exper- 
imental data in Fig. 2(a). Comparing these plots re- 
veals that, for an incident field of 24 kV/cm, a mobility 
of 4200 cm 2 /Vs (close to the reported value [24]) yields 
good agreement. If the majority of conduction electrons 
were scattered to the L valley, a mobility of 400 cm 2 /Vs 
would result. The agreement between experiment and 
simulation suggests that at 160 kV/cm this a reasonable 
approximation. 

Above 160 kV/cm, IMI becomes increasingly impor- 
tant. This is because the SRR resonance is strength- 
ened from the IVS-induced conductivity decrease. At 
160 kV/cm, simulations indicate in-gap electric field en- 
hancement by a factor of ~15 on resonance, correspond- 
ing to an overall peak electric field enhancement of ~ 
4 in the time domain. This is in good agreement with 
previously reported values [TS]. Thus, with an incident 
field strength of 160 kV/cm, the in-gap field is estimated 
to be ~0.64 MV/cm, which is sufficient for IMI in which 
THz-accelerated conduction band electrons collide with 
valance band electrons to result in carrier generation (in 
the L valley) [T5J [T3] - As a simple estimate, the motion 
of electrons under an external field is given as 

d(v(t)) = eE(t) (v(t)) 

dt m e t(E) 1 ' 

where (v(t)) is the average electron velocity, E{t) is the 
THz field, m e is the effective electron mass, and r is the 
relaxation time [25] , For n-type GaAs, IMI is favorable 
for electrons in the L valley when the kinetic energy of 
an electron reaches a threshold energy, 8th = 2.1 eV [2"6] . 
Model calculations show that with an initial carrier den- 
sity of 1 x 10 16 /cm 3 , an incident field of 400 kV/cm (1.6 
MV/cm in gap field) can increase the carrier density up to 
~ 10 17 /cm 3 . This process leads to an in-gap conductiv- 
ity of ~ 6 (ficm) -1 . The calculated conductivity arising 
from IMI leads to reasonable agreement between experi- 
ment and simulation (Fig. |2jb) and (d)). For agreement 
between experiment and simulation at the highest fields, 
an average in-gap electron mobility of 240 cm 2 /Vs was 
used, consistent with increased carrier-carrier scattering 
at higher densities and in agreement with the lower end 
of the L valley mobility previously reported in optical 
carrier generation studies [20l [21] . 

In contrast to doped GaAs, semi-insulating GaAs (SI- 
GaAs, carrier density about 1 x 10 7 /cm 3 ) with its low 
conductivity of 1 x 10~ 8 (ilcm)^ 1 yields a distinct SRR 
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FIG. 3: (a) Experimental electric field transmission as a func- 
tion of frequency for various incident peak electric fields for 
SRRs on SI-GaAs. (b) Simulated transmission for SRRs on 
SI-GaAs as a function of frequency for various in-gap con- 
ductivities. The dashed curve is the predicted response if the 
entire lateral area of the GaAs substrate (to a depth of 0.5 
fim) has a conductivity of 20 (f2cm) _ . This clearly disagrees 
with experiment suggesting that the carrier multiplication oc- 
curs only in the capacitive gaps of the SRRs. The inset shows 
the simulated on-resonance field enhancement. 



resonance and THz field enhancement even at low inci- 
dent field levels. The inset of Fig. [3] shows the SRR 
unit cell used to fabricate arrays on SI-GaAs. The size 
and period of the SRR resonators are the same as the 
resonators on the doped GaAs, thereby maintaining the 
resonance roughly at the same frequency. A smaller SRR 
gap (1.3 /im) was used to increase the field enhancement 
given the difficulty in initiating IMI on SI-GaAs due to 
the small initial carrier density. At our free-space THz 
field strengths, IMI is ineffective as a carrier generation 
mechanism in SI-GaAs. However, with the additional 
field enhancement available from SRRs, we observed a 
nonlinear response on SI-GaAs that we attribute mainly 
to IMI. 

Figure [jjja) shows the experimentally measured trans- 
mission of metamaterials as a function of frequency for 
various values of Ej„. At the lowest field (100 kV/cm), 
the response is still in the linear regime with a charac- 
teristic dip in the transmission at the LC resonance fre- 
quency. As the incident field increases, the metamate- 
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rial transmission increases and broadens with negligible 
change at lower frequencies. We note that even at the 
highest incident fields (Ej n = 400 kV/cm), no nonlin- 
ear transmission changes were detected on bare SI-GaAs. 
Similar to the nonlinear high field results (> 160 kV/cm) 
on n-type GaAs, the results of Fig. 3(a) suggest a large 
increase in the in-gap conductivity, indicating a carrier 
density increase via EVIL EM simulations were performed 
by changing the conductivity of the GaAs in the hor- 
izontal gaps as shown in the inset of Fig. 3(b). The 
plot shows the transmission change upon increasing the 
in-gap conductivity from 10 -8 (ficm) -1 to 20 (Ocm) -1 . 
The simulations show a trend that is markedly similar to 
the experimental results. 

The inset of Fig. [3^b) shows the field within the gaps 
is enhanced by over 20 times on resonance correspond- 
ing to a peak time-domain field enhancement factor of 
7. Impact ionization calculations show that with an ini- 
tial carrier density of 1 x 10 7 /cm 3 , an in-gap field of 2.8 
MV/cm (corresponding to our largest incident field, E in 
= 400 kV/cm) can increase the carrier density by more 
than ten orders of magnitude up to ~ 10 18 /cm 3 . Al- 
though we attribute the carrier multiplication primarily 
to IMI, Zener tunneling cannot, a priori, be discounted. 
Tunneling calculations |27j show that with an in-gap peak 
field of 2.8 MV/cm, Zener tunneling can generate a free 
carrier density of ~ 10 17 /cm 3 which is one order of mag- 
nitude smaller than that reached by IMI. This indicates 
that carrier generation may involve both processes. 



*3f 6 

tr 



SRRs on n-type GaAs 



SRRs on SI-GaAs 



(a) jfl 


24 


72 




100 




160 


Ej n (kWcm) ^^^l 




230 




— 300 K 




350 




385 Vj 




— 400 \/ 





-H.20 
DC .„ 



(c) / 




Ej n (kV/cm) 




100 




234 


300 




353 




388 
— 400 






50 150 250 350 
Incident E-field (kV/cm) 



100 200 300 400 

Incident E-field (kV/cm) 



FIG. 4: (a) Effective real permittivity e re (w) for SRRs on 
doped GaAs for various incident electric fields, (b) e re (uj) as 
a function of incident peak electric field strength at 0.82 THz 
for SRRs on doped GaAs. (c) e re {ui) for SRRs on SI-GaAs 
for various incident electric fields, (d) t r ^{oS) as a function of 
incident peak electric field strength at 0.85 THz for SRRs on 
SI-GaAs. 



One consequence of embedding metallic resonators into 
nonlinear dielectrics is the enhancement of the effective 
nonlinearity by 2 to 4 orders of magnitude (in compari- 
son to the bare nonlinear material) because of the in-gap 
field enhancement [35]. For our metamaterials, the non- 
linearity results in a change of the oscillator strength of 
the LC resonance. Because of dispersion, this leads to 
changes in the effective permittivity and refractive index 
at frequencies away from resonance. The effective real 
permittivity (e re ) is shown as a function of frequency for 
various incident electric fields for SRRs on n-type GaAs 
(Fig. ga)) and SI-GaAs (Fig. Qc)). These plots were 
obtained by performing parameter extraction (for a cubic 
unit cell) using the full-wave simulations [29H31] of Fig. 
gc), (d) andgb) and subsequently mapping the in-gap 
conductivity to the corresponding incident peak electric 
field required to achieve that conductivity. It is clear that 
the field-dependence of e re is different for the two cases. 
For the SRRs on n-type GaAs, the off-resonance e re at 
0.82 THz changes from positive to negative owing to IVS 
in the doped thin film. Above 160 kV/cm, e re increases 
from negative to positive values due to the onset of IMI 
(Fig. 4(b)). For SRRs on SI-GaAs, IMI increases e re 
monotonically from negative to positive values at 0.85 
THz (Fig. 4(d)). 

We have demonstrated that metamaterials enable 
strongly enhanced nonlinearities at terahertz frequen- 
cies. The present study illustrates strong THz-induced 
changes in the carrier mobility and density, yielding con- 
ductivity changes up to ten orders of magnitude and large 
changes in THz transmission. As our proof-of-principle 
results demonstrate, it will be possible to create a host 
of nonlinear THz metamaterials which includes, as but 
one example, nonlinear absorbers for saturable absorber 
or optical limiting applications. 
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